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Many one-dimensional quantum systems, in particular interacting electron and spin systems, can 
be described a Luttinger liquids. Here, some basic ideas of this picture of one-dimensional systems 
are briefly reviewed. I then discuss the effect of interchain coupling for a finite number of parallel 
chains. In the case of spin chains coupled by exchange interactions, the low-energy properties are 
radically different according to whether the number of coupled chains is even or odd: even number of 
chains have a gap in the spin excitations, whereas odd numbers of chains are gapless. The effect of 
interchain tunneling is analyzed for two and three coupled chains of itinerant fermions: for repulsive 
interactions, the two-chain system is "universally" found to be a d-wave superconductor, with a gap 
£^ in the spin excitation spectrum. On the other hand, for three chains the ground state depends both 
on the boundary conditions in the transverse direction and on the strength of the interactions. Weak 
repulsive interactions in all cases lead to dominant superconducting pairing of d-type. An example 
of a three-leg spin ladder with a spin gap is proposed. A general scheme to keep track of fermion 
anticommutation in the bosonization technique is developed. 
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<■> ■ I. INTRODUCTION 

> ' 

^ . Much of our current understanding of the physics of interacting electron systems, in partic- 

ular in metals, is based on Landau's Fermi liquid theory: |1| the properties of the interacting 
system (the Fermi liquid) closely resemble those of a system of noninteracting fermions (the 
Fermi gas). Interactions only lead to finite renormalizations of physical properties (specific 
heat, spin susceptibility, ...), as well as possibly to symmetry breaking phase transitions. Fun- 
damentally, the properties of a Fermi liquid are all related to the existence of quasiparticle 
excitations whose lifetime diverges as 1/uj 2 when the excitation energy uj goes to zero. This 
leads to universal behavior of correlation functions at low energies. 

A radically different situation is found in one-dimensional interacting systems: due to the 
restricted motion along only one direction in space (or equivalently, the not very restrictive 
conditions on scattering phase space) perturbations easily propagate coherently. This leads in 
particular to the fact that all the low-energy excitations are collective modes, namely long- 
wavelength fluctuations of the density (or spin density), and no quasiparticle type elementary 
excitations exist. Correlation functions then have nonuniversal, interaction-dependent power- 
law behavior at low energies, and for the case of spin-1/2 fermions the interesting phenomenon 
of spin-charge separation arises. 

In the following section I will review some important results from the theory of spinless 
interacting fermions in one dimension and their application in particular to spin-chain systems. 
Both isolated or coupled parallel chains will be considered. Most of these results are well-known, 
but they will be useful for the understanding of section JTT, where the case of spin-1/2 fermions 



will be considered and recent results on the coupled chain problem will be reported. For a 
general description of the subject, the reader is referred to review articles. 



II. THE SPINLESS LUTTINGER LIQUID: SPIN CHAINS 

A. Luttinger model and Luttinger liquid 

The basic ingredient of Luttinger liquid theory is the existence of an exactly solvable model 
of interacting fermions in one dimension, the Luttinger model. [pJU In this model the kinetic 
energy is essentially that of relativistic one-dimensional fermions, described by the Hamiltonian 

H = v F J2i( k - k F )a\ k a + k + (-k - k F )a ] _ )k a_ k } , (1) 

k 

where I use a + -(a„-)operators for right- (left-) moving fermions, and Vp is the Fermi velocity. 
Here the /c-summation runs from -oo to oo. In the noninteracting ground state all negative 
energy states are filled. Note that kp can be removed by a simple unitary transformation. 
Interactions between fermions are of the forward scattering type, (kp; —kp) — > (kp; —kp) or 
(kp; kp) — > (kp; k F ). The Hamiltonian describing these processes takes the form 

H im = ^E{M?K(?)p-(-?) + 94(q)[p+(q)p+(-q) + p-(-q)p-(q)]} ■ (2) 

Here, (72(9) and g^(q) are the Fourier transforms of a real space interaction potential, p±(q) 
are the Fourier components of the fermion density operator for right-(+) and left-(-)-going 
fermions, and L is the length of the system. 

The model defined by eqs.(HD and @j can be solved exactly by the bosonization technique. || 
One introduces a scalar boson field <p(x) and its conjugate momentum density H(x), satisfying 
canonical boson commutation relations: 

[<t>(x),IL(y)] = i6(x-y) . (3) 

4>(x) is related to the local fermion density via 

d x (j)(x) = -Ti[p(x) - p ] = -n[tp+(x)i) + (x) + ipl(x)ip_(x) - p ] . (4) 

Here po is the average particle density, and the ip± are standard fermion field operators. In 
terms of 6 and n the full Hamiltonian becomes 



H = H + H int = / dx 



(5) 



This is obviously just the Hamiltonian of an elastic string, with the eigenmodes corresponding 
to the collective density fluctuations of the fermion liquid. It is important to notice that these 
collective modes are the only (low-energy) excited states, and that in particular there are no 
well-defined single particle excitations. The parameters in (El) are given by 



u=[(vp + g 4 /(2n)) 2 -g 2 2 /(2nf]^ , K 



2-KVp + g 4 - g 2 
2nv F + c/4 + g 2 



(6) 



The energies of the eigenstates are u(q) = u\q\, i.e. the parameter u determines the phase 
velocity of the low-lying excitations. 



The final ingredient for the following is the expression for the single-fermion field operators 
01 
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ip±(x) = —;= a± k e lkx = lim —j=rj± exp [±ik F x =F i(j>(x) + i9{x)\ . (7) 



Here d x 6(x) = 7rn(x), a is a short-range cutoff parameter, and the r]± are Majorana ("real") 
fermion operators introduced to guarantee proper anticommutation between ip + and ■?/>_. || 
They satisfy the anticommutation relation 

[Vr,Vs\+ = 2<5 r , s , (8) 

implying in particular (r/ r ) 2 = 1. They can be represented by standard (Dirac) fermion operators 
c r as r) r — c\ + c r . Note that there is just one isolated fermionic degree of freedom per branch, 
and that these degrees of freedom do not appear in the bosonized Hamiltonian. 

Using eq.(^) correlation functions can now be easily calculated. One has the general result 
for ground state expectation values of the Hamiltonian (^) 

T](m,n) = ^(m 2 K + n 2 /K) . (9) 

Here x{ x i T ) ls the angle between the vector (x, ur) and the x-axis, and r is the Matsubara imag- 
inary time. Thus, quite generally, correlation function decay asymptotically with nonuniversal 
power laws which exponents depending on the strength of the interaction via the coefficient 
K, eq.(|6|). After Fourier transformation this of course gives nonuniversal power laws at low 
energies and long wavelengths. As an example consider the pairing correlation function 

(V4(x,t)V>1(z,t)V^0,0)V^(0,0)) « (e 2ie ^) e - 2i9(0 ' 0) ) oc (x 2 + u 2 t 2 )~ 1/k . (10) 

As expected in a one-dimensional model, this correlation function never shows long-range or- 
der, however, for attractive interaction (corresponding to K > 1) where standard BCS theory 
would predict long-range order, the decay is slower than in the noninteracting case, indicat- 
ing enhanced pairing correlations. Similar results can be found for a large variety of other 
correlation functions. 

The above is exact for the Luttinger model. On the other hand, one would expect that a 
model given by eqs.(|l|) and (^) should also be appropriate for the low-energy properties of, for 
example, models for interacting fermions moving on a lattice or with nonlinear energy disper- 
sion. This is certainly correct in the weak-interaction limit where all the low-energy physics 
is dominated by processes near the Fermi level and the structure of states far from the Fermi 
level is largely irrelevant. By continuity as long as no phase transition occurs as a function 
of increasing interaction strength, eq.(^) then also is the effective low-energy Hamiltonian for 
stronger interactions. The coefficients u and K of course depend in a more or less complicated 
fashion on the interaction, but the basic fact that all the low-energy physics is determined 
entirely by u and K remains. In particular, scaling relations between exponents of different 
correlation functions are always the same. This reasoning can be made quantitative using renor- 
malization group arguments to show that (almost) all the effects neglected (band curvature, 
finite cutoffs, momentum dependence of the interaction,...) are irrelevant and therefore only 
can lead to quantitative renormalizations of u and K. JTOj] Physical systems to which this kind 



of description applies include, apart from interacting electrons, also one-dimensional spin and 
interacting boson systems, as well a the edge states of quantum Hall devices. These system, 
where the Hamiltonian is correct for low-energy properties but in general not for high-energies 



(of order of the bandwidth or higher), are now referred to as Luttinger liquids. |T(J In all cases, 
the coefficients u and K, sometimes appropriately generalized to include the effect of electron 
spin, determine low-energy correlation functions, and therefore transport properties, the effect 
of extra perturbations like random potentials, and low energy spectroscopic properties. 



B. Spin chain problems 

A particularly interesting example of the use of the Luttinger liquid concept is provided by 
the spin-1/2 antiferromagnetic chain, described by the spin Hamiltonian 



H = ^(SfSf +1 + S?SV +1 + ASfSf +1 ) 



(11) 



Here Si = (Sf, Sf, Sf) is a spin-1/2 operator acting on site i, A is an anisotropy parameter that 
allows one to treat the antiferromagnetic (A = 1), the ferromagnetic (A = —1), and general 
anisotropic cases. The model has been solved exactly by Bethe's celebrated "ansatz", |TT| , |r2"[ 
but many of its properties have remained difficult to describe until much later. 

The Luttinger liquid concept in fact is very useful for describing asymptotics of correlation 
functions for this model. One starts by the Jordan-Wigner transformation JT3|] from spin to 
fermion operators (a i; a|): 



5+ 



Sf + iSf = a J exp 



2-1 

in '^2 a ] a j 

.7 = 1 



Sa 



1 



(12) 



In the fermionic language, the SfS~ +l interaction leads to a nearest neighbor hopping term, 
and SfSf +1 leads to a nearest neighbor interaction. Following the logic outlined above, one 
then finds eq.([5|) as the effective low-energy Hamiltonian, with in particular K = 1 for A = 0. 
The dominant (i.e. most slowly decaying) contributions to the spin operators are 
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-l) i cos(2(/)(m)) 



(13) 



where a is the lattice constant. From eq.(||) the longitudinal and transverse spin-spin correlation 
functions then decay as x~ 2K and x" 1 '^', respectively. This immediately implies that for 
the isotropic antiferromagnet (A = 1), where spin rotation invariance imposes equality of 
longitudinal and transverse correlations, one has 



(Sj(t) • S ) 



■ ln^V + Mj 



x 2 + (ut) 



x 



ja 



(14) 



Here the logarithmic factor p4]-4T6| comes from a marginally irrelevant umklapp operator |T 
(oc cos(40) in bosonic language) which has been neglected up to here. We thus have K = 1/2 
for the isotropic antiferromagnet. Eq.flTJI) and its generalization to finite temperature can be 
used to obtain an analytic expression for the magnetic neutron scattering cross section [18 
which is found to be in excellent agreement with experiment. |i~9|j20[| 

For general values of A (obeying |A| < 1) K can be determined noting that in eq.(|5D the 
coefficient u/K determines the inverse susceptibility, which can be obtained from the exact 
Bethe Ansatz solution of the spin model. One then obtains [pLO 



K = 



71 

2 arccos(— A) 



vrVl - A 2 
2 arccos A 



(15) 



C. Coupled spin chains 



Investigating models of coupled parallel chains is of interest for a number of reasons: (i) 
guasi-one-dimensional antiferromagnets always have some form of interchain coupling, usually 
leading to three-dimensional ordering at sufficiently low temperatures; |22|-p4| (ii) there is a 
number of "spin-ladder" compounds containing a small number of coupled chains; p5| , |26[| (iii) 
coupled spin-1/2 chain models can be used to describe higher spin quantum numbers. [p7| 

Consider N coupled spin-1/2 chains with spin degrees of freedom Sj, j = 1..N, described 
by the Hamiltonian 

N 

H = ^2H(Sj) + ^2\ jk H c (Sj, S k ) , H c (Sj, S k ) = ^2Sj ti ■ S k ^ , (16) 

j=l j<k i 

where i labels sites along the chains, j,k label the chains, and H(Sj) is of the form (|TT|). For 
Xj k = — 1 the ground state of this model is exactly that of the spin-iV/2 chain, each site being in 
a state of total spin N/2, |27[] but the model can be considered for general A, both ferromagnetic 
(A < 0) and antiferromagnetic (A > 0). Performing now the Jordan- Wigner transformation for 



the Sj separately and going to the boson representation the Hamiltonian becomes [|18 
H = I dx 



+ i 1 , 2 [ dx iKjk cos(2(0j + 4> k )) + \ 2 ,j k cos(2(0j -<f>k)) + A 3 jfc cos(0j - 6 k )} (17) 
ytoi) j<k J 

Here <fi = (<f>i, 2 , .., (f>jy), 4> = J2j4>j/VN, the constants u,K,K all depend on the different 
constants in the original Hamiltonian, and the Xx^-jk are au proportional to Xj k . 

Elementary power counting, using the result eq.(|l4]) for the spin correlations, shows that 
the coupling term H c always is a strongly relevant perturbation. Consequently, an explicit 



renormalization group calculation |18[ shows that either A2 or A3 always scale to strong coupling, 
i.e. the "relative" degrees of freedom <f>j — (fi k all acquire a gap. In particular, for not too strong 
anisotropy, the A3 operator dominates, giving rise to long-range order in the 6j — 9 k and 
correspondingly exponential decay of the <pj — (p k correlations. Integrating out these massive 
degrees of freedom an effective Hamiltonian for the "global" mode is found: 

wKn(x) 2 + ^(dJ) 2 + gcos(iiVN^-—hd x ^\ , (18) 

TlK 7T J 

where [i = 2 for even N and n = 4 for odd N, the coefficients u, K, g are renormalized, and h 
is an external magnetic field applied along the z direction. Similarly, the leading contribution 
to spin correlations comes from the operators 

S + {x) oc e i7rx e~ i§/VN , S z {x) oc e inx cos(2v / iV0) , (19) 
where the second equation applies to odd only. 




From this a number of important conclusions can be drawn. [|18| We first notice that massless 
excitations and the corresponding slow algebraic decay of correlation functions are only possible 
if the cos term in eq.(|18D is irrelevant. Moreover from eqs.(^) and (^) it follows that spin 
correlation functions are isotropic only if K — 1/{2N), implying a decay as 1/x, as in the 
5 = 1/2 case. For the case of odd N (equivalently, for ferromagnetic A, for half-odd-integer S) 
this is indeed the correct behavior: for fi = 4 the cos term is marginally irrelevant. Thus both 
antiferromagnetic spin chains with half-odd-integer S and for odd numbers of coupled 5 = 1/2 



chains massless behavior is predicted, with correlations asymptotically decaying like those of a 
spin-1/2 chain. |18|] There is both numerical [28-|31Jl and experimental [^] evidence that this 
is correct. 

On the other hand, for even N ( equivalent ly, integer S) the cos term in eq. (|T8|) is strongly 
relevant and therefore generates a gap A s in the spin excitations. For the integer-S 1 spin chains 
this is the well-known and verified Haldane prediction f32fl , but there also is a gap for any 
even number of coupled chains. |fIB| , pB[ The gap implies exponential decay of spin correlations. 
Numerical [|30|j3TJ| and, at least for N = 2, experimental work [^5j confirms this picture. Another 
prediction, again valid both for integer-S* antiferromagnets and even numbers of coupled chains, 
concerns the effect of an applied magnetic field: p^j3^j35| as long as the field is smaller than 
a critical field h c oc A s , the ground state is unchanged and has zero magnetization. However, 
beyond h c the magnetization is expected to increase as M oc \fh — h c . Experiments on S — 1 
antiferromagnetic chains confirm this prediction. ]36| 

A natural question left open by the above considerations concerns quasi-one-dimensional 
antiferromagnets like KCuF 3 |22|,|23| or Sr 2 Cu03, fl24]j37|l which can be considered as the N — > oo 
limit of the above model. One clearly expects (and observes) true antiferromagnetic order at 
sufficiently low temperatures, at first sight in contradiction both with the exponential decay 
of spin correlation predicted for even N and the universal 1/x law for odd N. However, one 
should note that on the one hand the correlation length of the even-iV systems is expected 
to increase quickly with increasing N, and that on the other hand the 1/x correlation law 
of the odd- A systems also is expected to be only valid beyond a correlation length £(N) 
which increase with N. In the thermodynamic limit N — > oo this then is perfectly consistent 
with the existence of long-range order. Theoretical treatments of magnetic order in quasi- 
one-dimensional antiferromagnets can be based on a mean-field treatment of the interchain 
interaction P5|,|5P| which gives quantitative predictions for systems like KCuF 3 or Sr 2 Cu0 3 . 



III. LUTTINGER LIQUID WITH SPIN: (QUASI ) ONE DIMENSIONAL 

CONDUCTORS 

A. One chain... 

The generalization to fermions with spin of the method of the previous section is relatively 
straightforward: |2|] one introduces one 0-field for each spin orientation, and then uses the 
charge-(p) and spin-(a) fields (j) Pj(T = (0-j- ± 0|)/v2- The dynamics of theses fields is governed 
by separate charge- and spin-Hamiltonians, each of the form (^), with constants u Pt(J , K p>a . 
In addition, backward scattering leads to the existence of an extra term g\ cos(v / 80 CT )- For 
repulsive interactions (gi > 0) this term renormalizes to zero, and I will concentrate on this case 
in the following. By arguments similar to those used for the spin chain problem, spin rotation 
invariance requires K a = 1, and then there are three interaction dependent constants u p , u a , K p 
left in the effective low-energy Hamiltonian. These can be determined either perturbatively, 
giving rise to expressions like eq.(^), or can be determined from exact solutions and numerical 



finite-size diagonalization. [^0| , ^T| For a noninteracting system u p = u a = vp, K p = 1. 

A remarkable feature is the separation of spin- and charge degrees of freedom: || the Hamil- 
tonian is a sum of spin- and charge part, and, from eq.(^) operators and therefore correlation 
functions become products of a spin- and a charge part. As spin and charge excitations prop- 
agate with different velocities, this gives rise to the interesting phenomenon of spin-charge 
separation, leading for example to a complicated double-peak structure in the single-particle 
spectral density 'Q,^] instead of the single quasiparticle peak of a Fermi liquid. Charge and 



spin-modes propagating at different velocities have indeed been observed in experiments on 
quantum wires. [44,45] Of course, as in the spinless case, no broken-symmetry ground states 



occur, but fluctuations of various types of order are enhanced. In particular, for purely repulsive 
interactions, antiferromagnetic (equivalently, spin density wave) fluctuations dominate. 



B. ...two chains... 

It is clearly of interest to see what happens to the peculiar one-dimensional behavior when 
one puts chains in parallel. This question is of relevance for the understanding of quasi- 
one-dimensional conductors, doped spin ladders, few-channel quantum wires, and generally 
for the understanding of possible non-Fermi-liquid behavior and correlation-induced super- 
conductivity in higher-dimensional solids. Of particular interest is the effect of an interchain 
single-particle tunneling term of the form 



Hij 



t± / dx(ipl si ifj rsj + h.c.) 



(20) 



where ip rs j is the fermion field operator for right (r = +) or left (r = — ) going particles 
of spin s on chain i. Simple scaling arguments |4B| lead to the "phase diagram" shown in 
fig.0. The dashed line represents the crossover below which single-particle tunneling becomes 
strongly relevant and below which one thus expects Fermi liquid like behavior (an alternative 
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FIG. 1. Qualitative phase diagram in the temper atwce-K p plane for Luttinger liquids coupled by 
interchain hopping. K± = 3 ± \/8- 



interpretation is due to Anderson [47,[4"5| ). The full lines indicate where two-particle or particle- 
hole hopping becomes relevant. The most plausible interpretation is that this is the temperature 
where three-dimensional long-range order of some type sets in. A more detailed discussion of 
this in quasi-one-dimensional systems (a thermodynamically large number of chains) has been 
given in this conference. |49|j50[1 

As a model for potential doped spin-ladder systems, as well as a first step towards a many- 
chain system, one can study the two-chain case. |5l]-^| The tunneling term, eq.(^Cj), then leads 
to a splitting of the single-particle bands into symmetric and antisymmetric combinations which 
I label by transverse wavenumbers k± = 0, tt. Now each k± mode can be bosonized separately. 
Introducing the linear combinations <p u ± = ((j) u o±(j) un ) / \/2 {y = p, a) the Hamiltonian (including 
t±) then takes the form 



H = H + H int i + H, 
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(21) 



fl'oooo + Tfl'o^lrOj an d S'ifecd is the coupling constant for an interaction scattering two 



Here 

particles from fcj_-states (a, b) into (d, c). The signs of the different interaction terms in eq.(21) 
have been determined following the reasoning explained in the appendix. Here I consider a 
case where there is only intrachain interaction, implying that all bare coupling constants are 
independent of k±. 

For the pure forward scattering model (g± = 0) the only nonlinear interaction (^qoItt) scales 
to infinity, leading to a gap in the (p— ) modes and in half of the (a— ) modes. The remaining 
(er— ) modes are protected by the duality symmetry under <p a _ <-» 9 a _, the (er— ) sector is 
in fact a critical point of the Ising type for g\ = 0. [53] For repulsive interactions here the 
dominant fluctuations are of CDW type, with decay proportional to r - < ^+ 2K p)/ 4 ( an d K 2 p = 
(tivf — gi + g\/2) / (tivf + g% — fl'i/2). This state can be label by he number of massless modes 
as ClSl^, where quite generally CnSm denotes a state with n massless charge and m massless 
spin modes. 



For nonzero g\ all interactions scale to strong coupling, only the total charge mode remains 
massless (C1S0), reflecting the translational invariance of the system, and all spin excitations 
have a gap. The physics in this regime can be determined looking for the semiclassical minima 
of the different cos terms in eq. fl2T|) . One then finds that the CDW correlations now decay 
exponentially, and for the interesting case g\ > 0, gi > gi/2, corresponding to purely repulsive 
interaction, the strongest fluctuations are now of "d-type" superconducting pairing, |5T| , |5"3"| 
with decay as r~ 1 ^ 2Kp \ Labeling this state as "d-type" seems appropriate because the pairing 
amplitudes at kj_ = and tc intervene with opposite sign. In real space, this corresponds to 
pairs formed by two fermions on the two different chains. Note that even for weak interactions 
where K p — > 1 this decay is very slow. The 4k p component of the density correlations also 
has a power law decay, however with an exponent 2K P , much bigger than the SCd exponent. 
5B],[5^] The full phase diagram in the g\— g 2 plane is shown in fig.0. For g\ < the diagram 



is identical to the single-chain case, however, for g± > the behavior is changed dramatically, 
and in particular superconductivity is predicted for repulsive interactions, for example for the 
Hubbard model which would be represented as g\ = g 2 in the present language. 

Remarkably, results basically identical to this weak coupling analysis can be obtained as- 
suming strong repulsive interactions in the individual chains, so that one is for example in 
the regime where interchain electron-hole pair hopping is more relevant than single particle 
tunneling, K p < 1/3 in fig.|TJ) . It is then more appropriate to bosonize the degrees of freedom of 
individual chains, rather than working in k± -space. Under renormalization one then generates 
an interchain interaction of the form [I6,|55| 
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where j, k now are chain indices. Remarkably, this term leads to properties identical to those 



found in weak coupling. [53| First, for #i = one has V = J±. Then Hj k is invariant under the 
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FIG. 2. Phase diagram of the two-chain model. The different dominating fluctuations are: SCd: 
"d-type" pairing; OAF: orbital antiferromagnetism; SS: "s-type" pairing; CDW: charge density wave. 
The critical lines g\ = and g\ = 2g 2 are of Ising type. 



duality <p a j <-> 9 a j, and one has an Ising critical theory. Secondly, for #1 > one finds V ^ J±. 
Then the fields (f> p i — <p p2 , 4> a i + 4>a2, and Q a \ — 9 a2 become long-range ordered and one is in 
a C1S0 state. Power law correlations again exist for SCd and density fluctuations, with 
the same scaling relation between the two exponents as in weak coupling. However, because 
now K p < 1/3 the Akp density fluctuations actually dominate. The equivalent results in the 
weak and strong coupling regime very strongly suggest that the two-chain model is in the same 
phase for weak and strong repulsion. This point is further supported by considering the "t-J 
ladder model" for strong interchain exchange. |56fl 

Numerical work on the two-chain model is in agreement with the existence of d-type pairing, 
|57|-|6Tf the evidence for the 4fcp density fluctuations is however inconclusive. Concerning 
experimental observation, one should notice that the SCd state becomes localized by weak 
disorder. 



C. ...three chains 

Given the differences between even and odd numbers of parallel spin chains discussed in 
sec.jry, it is clearly of interest to see what happens when one passes from two to three (and 
eventually more) chains. For the three-chain case it is important to consider boundary condi- 
tions in the transverse direction carefully: periodic boundary conditions are frustrating, and I 
therefore consider both periodic and open boundary conditions. The corresponding geometries 
and electronic structures near the Fermi energy are shown in fig.0. 

In the periodic case, there are 12 different coupling constants implying particles at the Fermi 
energy. These interactions do not satisfy the criteria for bosonization developed in the appendix. 
However, a one-loop renormalization group calculation shows that certain of these interactions 
scale to strong coupling, whereas others remain finite. It turns out that the interactions that 
do scale to strong coupling do also satisfy the criteria of the appendix, i.e. following the 
same logic as in the two-chain case the strong coupling regime can be understood based on a 
bosonized picture. Specifically, I then find: (i) if only forward scattering is concerned, one is at 
an Ising critical point similar to the two-chain case, in a C1S2\ state. Both CDW and SDW 
correlations decay exponentially; (ii) upon introducing backward scattering, all the massless 
spin modes disappear (CISD); (iii) for gi > and gi < 2g 2 the dominant fluctuations are of 
"cfci/'-type, where by this I mean that the pairing amplitude has opposite signs ant k± = ±27r/3 
and vanishes at k± = 0; for g\ > and g\ > 2g 2 the pairing amplitude has the same sign at all 
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FIG. 3. Structure of three parallel chains with periodic (top) and open (bottom) boundary condi- 
tions in the transverse directions. Full and dashed lines are the perpendicular and parallel hopping 
matrix elements, respectively. To the right the bands crossing the Fermi energy are shown schemat- 
ically. In the periodic case, the bands with k± = ±27r/3 are degenerate. In the open case bands are 
labeled by their transverse quantum number, as explained in the text. 



k±. In both cases, the decay exponent is l/(3K p ), with K p given by an expression analogous 
to the two-chain case. We note that at least for two and three chains, the exponent takes the 
form l/(NK p ), with N the number of chains. For arbitrary N this form can be obtained from 
a microscopically derived Ginzburg-Landau description. 

For open boundary conditions in the transverse directions, I label eigenstates on the three 
chains as |— ), |0), |+), as in fig.|3]. The eigenstates of the transverse hopping operators then are 
= (|+) - \-))/V2, |0> = (|+) + |0) + |-))/V3, and | + 1) = (|+) - 2|0) + \-))/VB. We 
now have 18 independent coupling constants. However, the renormalization group equations 
show that the antisymmetric | — 1) states nearly decouple from the |0) and |+1) modes, i.e. the 
problem basically reduces to that of uncoupled single- and double-chain systems. For g 1 = 
one then is in a C2S2^ state with Ising criticality, whereas for gi > one has a C2S1 state. 
The dominant correlations are a form of interchain pairing involving mainly particles in the 
|0) and | + 1) states. The decay exponent is l/(3K p0 ) + V(6-Kpi), approximately 1/2 for weak 
interaction Here K pQ refers to the total charge mode and K p i to the other massless charge mode. 
There are also subdominant SDW fluctuations coming from the |— 1) transverse states, with 
exponent 1 + l/(3K p0 ) + 2/(3K p i), close to 2 for weak interaction. The number of massless 
modes for g% ^ have also been determined by Arrigoni. f64jl 

The above is based on an analysis of the perturbative renormalization group equations 
and therefore valid for weak interactions. For strong intrachain interactions one can proceed 
similarly to the two-chain case: for K p < 1/3 the dominant interchain interaction, generated 
from the initial single-particle hopping, is of the form of eq . fl22|) . As the interaction constants 
renormalize to strong coupling, the semiclassical solution minimizing this term is <p P j — <p p k ~ 0. 
After the unitary transformation <pj = (j) a j/y/2, 8j = V29 a j Hjk then takes the same form as 
the interchain coupling term of the coupled spin chain model, eq. (|PT|) , and the analysis of that 
Hamiltonian can be taken over. Thus, for an odd number N of chains (and in particular for 
N = 3) one has one massless spin mode, and in addition a massless charge mode representing 
in-phase oscillations of all chains, the state thus is C1S1. The dominant fluctuations are of 
SDW type which decay as r~ l ~ Kp ' N . On the other hand, for all even N spin excitations have 
a gap that decreases with increasing N, and the state is labeled C1S0. Similar conclusions can 
be reached analyzing t-J ladders in the limit of large interchain exchange. 

The analysis of the preceding paragraph applies as long as the interchain exchange interac- 
tion is unfrustrated, e.g. as far as iV = 3 is concerned to open boundary conditions. For the 
case of periodic boundary conditions in the transverse directions the interchain exchange still 
scales to strong coupling, however the resulting model has quite different behavior: as far as 



the spin degrees of freedom are concerned, one essentially recovers a Heisenberg model on the 
lattice shown in the top part of fig.|3], with exchange constants J± within each triangle and J 
between triangles. In the limit J± ^> J, towards which the renormalization group scales, this 
model can be further transformed: in the ground state, each triangle is expected to be in one 
of the spin-1/2 states 

iTi?> = -^(im>+^im>+^ 2 im>) , iTL> = -i=(im>+^ 2 iUT> + ^im>) , 

|Ii?) = i=(|tlI) + ^|in) + ^ 2 |!IT)) , |U> = ^(|U!> + ^ 2 im> + c|UT» , (23) 

where u = exp(27ri/3). States with total spin 3/2 on a triangle are by J± higher in energy and 
thus are neglected. In this subspace, the effective Hamiltonian becomes 

HeS = | E St ■ S l+1 [l + 2(rfTf +1 + t?t? +1 )} , (24) 



where i label triangles, S% is a spin-1/2 operator acting on the first indices of the states (|23|). 
and the r, are Pauli matrices acting on the second index in (|23|) . The Hamiltonian (|24]) can now 
be analyzed using the Jordan-Wigner transformation. In the fermionic language there then are 
both umklapp and backward scattering process that lead to gaps in all excitation, i.e. all the 
spin excitations are gapped. For the original three-chain problem then only the total charge 
mode is left, the state can be labeled as C1S0. Interestingly, this is the same situation as that 
found in the preceding weak coupling analysis. Note also that we have here a counterexample 
to the general discussion of sec. [ll C| , namely, a three-chain system with a (frustration-induced) 
spin gap. 

Concluding the analysis of the three chain problem we may notice that, similarly to the 
two-chain case, for weak coupling we find dominant d-wave (or interchain pairing) fluctuations. 
However, for stronger intrachain interactions at least in the case of nonfrustrating boundary 
conditions in the transverse directions, the behavior is quite different from weak coupling, 
showing in particular dominant SDW fluctuations. This is quite different from the "universal" 
properties of the two-chain case. The principal reason seems to lie in the different mechanisms 
for interchain coupling: for weak interaction single-particle tunneling dominates, whereas for 
stronger interaction (induced) particle-hole process are more important. 



IV. CONCLUSIONS 

In this paper I've discussed a number of results, mainly analytical, on the effect of differ- 
ent forms of interchain coupling on the Luttinger liquid behavior of strictly one-dimensional 
systems. As far as spin chains are concerned, the most spectacular result is the "oscillation" 
between even and odd numbers of chains, reminiscent of (and formally related) the Haldane 
phenomenon in spin-S antiferromagnetic chains. There is both experimental and numerical 
evidence for this behavior. For conducting chains, the situation seems to be more complicated: 
for very strong intrachain repulsion, similar "oscillations" as for the spin chains are predicted, 
however, the number of massless modes also depends on the strength of the interaction. At least 
for weak interaction, there are dominant d-type pairing fluctuations both for two and three 
chains. Numerical work on the two-chain system is mostly in agreement with the analytical 
results. These results in principle apply directly to doped spin ladders, on which experimental 
results however are rather scarce at the moment, and to few-channel quantum wires. In both 
cases are more detailed understanding of the effect of disorder is needed. 
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APPENDIX A: WHEN TO BOSONIZE IN PEACE 

Using eq. (0) and its generalization to cases with spin and other "internal" degrees of freedom 
like perpendicular momentum indices in coupled chain problems, a typical fermion interaction 
term becomes 

^a^aV^a = VaVpV-yVs x (boson operators) = h a/3y s x (boson operators) , (Al) 

where the second equality defines h al 3 y s- This operator, responsible for taking into account 
fermion anticommutation properly, is nevertheless carefully passed under the rug in the vast 
majority of the literature, thus leaving a purely bosonic Hamiltonian to be considered, as implied 
in the term "bosonization" . I will here investigate under which conditions this is allowed. 

The relevant situation is that all the indices a,(3,j6 in h a ^ are different from each other. 
Otherwise the anticommutation rule 

[r ]r , Vs ] + = 26 r ^^r ] 2 r = 1 , (A2) 
allows to simplify h a pjs- I will therefore only consider the general case. First note that 

h 2 a0jS = 1 , (A3) 
h af3l s thus has eigenvalues ±1. Secondly, 

[ha/3-,6, h K \^ v ]± = , (A4) 

where according to whether an even or odd number of pairs of indices taken from the two sets 
(01,(3,76) and (k, A, \i, v) are equal the commutator (even case) or anticommutator (odd case) 
is to be used. Finally, permutation of indices leads to sign changes: 

h a p~f5 — ~h/3a 7 s = —h ai ps — —h a /3s-/ (A5) 

It is now clear that if all the Ti's occuring in a given Hamiltonian commute, they can be 
simultaneously diagonalized, which means that it will be possible to replace each of the Ti's by 
±1, leading to a purely bosonic Hamiltonian. This clearly is the case if all the h's occuring have 
an even number of indices in common. In the opposite case some of the h's do not commute, 
therefore can not be simultaneously diagonalized and not be eliminated from the Hamiltonian. 
Bosonization then is not possible. 



As a simple example consider the single-chain Luttinger model with spin, sec.|HI[ the 
four allowed values of the discrete indices are 1 = (+,f),2 = (+,!), 3 = (— ,f),4 = (— , | 
). Consequently, only ^1234 can occur, and according to the eigenvalue chosen the backward 
scattering interaction takes the form ±g 1 cos(\^8(j) a ) . The choice of eigenvalue of /11234 affects 
however the expressions for correlation functions: for example ^1234 = ±1 implies 771773 = ±772774, 
and consequently the 2k p charge density operator contains either a factor cos(v / 20 CT ) (plus sign) 
or sin(v / 20 cr ). A similar discrete "gauge covariance" exists of course for all correlation functions. 

In more complicated cases like the two-chain problem, more then one /i-operator occurs. 
Even if they all commute, as is the case for the two-chain problem, additional constraints on 
the permissible eigenvalues of the Tt's exist due to the existence of relations of the type 

hap-ysh^x^h^p^j- = ±1 , (A6) 



and similar relations involving more than three h's. However, a discrete gauge freedom of the 
type mentioned above often remains. For the particular case of fermions with an internal SU(N) 
symmetry, |9[] bosonization can be performed without problem and all the /i-operators can be 
given eigenvalue +1, a fact not noticed in the original work. 
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